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ABSTRACT 

In this paper, we present a briefly overview of the protoplanetary 
disks in the Orion Nebula, incluiding some astrobiological aspects and 
an Ho Fabry-Perot study of 16 of them. We found that Fabry-Perot 
interferometry constitutes an effective technique for the detection 
of proplyds. We also report heliocentric systemic velocities for 
the proplyds 82-336, 158-323, 158-326, 159-350, 161-314, 161-324, 
163-317, 166-316, 167-317, 168-326, 170-337, 176-325, 177-341, 
180-331, 197-427 and 244-440. The velocities were measured between 
22-38 km s“^. 


Subject headings: CIRCUMSTELLAR MATTER — HII REGIONS — 
ISM: INDIVIDUAL (ORION NEBULA) — STARS: FORMATION — ISM: 
KINEMATICS 



1. INTRODUCTION 


The proplyds (PROtoPLanetarY DiskS) are photoevaporating circumstellar disks 
around young stellar objects (YSOs). They are composed of: 1.- a circumstellar disk around 
the YSOs, 2.- a head with a bright cusp, 3.- stand-off shocks and 4.- a tails. Their shapes 
are created when the UV photoionizing flow from a hot (O spectral type) external star 
interacts with a photoevaporated flow from the circumstellar disk: the photoionizing flux 
(2000 > A > 912 A) infringes on the disk producing a slow (~ 3 km s“A photodissociated 
flow. The flow passes through a weak D-type ionization front (IF) and is accelerated to 
transonic velocities (~ 10 km s“^) producing a weak shock that lies just inside the IF, 
generating a slight density and column density enhancement that is seen in silhouette 
(head). The electron density (ug) at the IF has values of ~ 10^-10^ cm“^. The greater values 
correspond to the inner regions. The 10® cm“^ density is reached just outside of the IF. The 
lie decreases with the distance from the YSO as r“^ due to the quasi-spherical geometry of 
the photoevaporated flow. The bright proplyd cusp corresponds to the photoevaporated 
flow from the head of the proplyd driven by the ionizing photons from the hot external star. 
Therefore, the head is composed of check ionization fronts and photoablated material from 
the circumstellar disk. The photoablated material expands outwards from this IF. When its 
ram pressure is balanced with the ram pressure of the photoionizing flow or stellar wind of 
the hot star, an arc-shaped stand-off shock forms, which is visible in Ha and in [OIII] lines. 
In these arcs, Ug^lO"^ cm“^. If the stellar infringing outflow is subsonic the Ha-1-[OIII] arcs 
are not formed. Finally, the tails are generated by the photoevaporated flow driven by the 
diffuse UV photon radiation from the nebula. These tails have Ug ~10® cm“® near to the 
head and ~10"^ cm“® at the tip. Recently, microjets have been observed in some proplyds 
[22, 3]. Detailed information about this scheme for proplyd formation is presented by [1], 
see also Figure 5 of [14]. More detailed, the actual proplyd formation schemes propose that 
the proplyds are formed and detected due to the photoevaporation of the protoplanetary 
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disk by the intense UV radiation of an external star (in this case 0^ Orionis C and, for 
some other proplyds, together with 0^ Orionis A). It is assumed that the protoplanetary 
disk is surrounded by a neutral envelope probably formed by a slow photodissociated wind 
from the disk. It is also assumed that an IF, forms in this envelope and that the newly 
ionized gas flows away from the proplyd with an initial velocity vq. Pressure gradients 
in the ionized gas accelerate the flow away from the IF. The amount of the acceleration 
depends on the type of the IF, being maximal when the flow leaving the IF surface is sonic 
(D-critical IF). The density in the ionized photoevaporated flow is highest at the point on 
the IF closest to the ionizing star, no, and decreases towards the sides. Assuming the flow 
to be isothermal and rg the radius of the IF to the point of maximum density, then the 
photoevaporated flow initial velocity, vq, is related to the disk’s photoevaporation mass-loss 
rate via the relation: M = dvr ro^ vq uq m/ (where m/ = 1.35 m^^) or [M/10“’^ yr~^] 

= 0.44 [ro/10^^ cm]^ [vo/km s“^] [no/cm“^] and the life time, r, of the protoplanetary disk 
is given by r ~ yidisk/M [14, 16]. Thus, to obtain the radial velocity prohles of proplyds 
is important because it allows us to constrain the photoevaporated flow velocity. This 
quantity gives us an estimate of the mass-loss rates of the protoplanetary disks and of the 
life time of those disks. The scheme described above is also presented by [13, 11, 14]. and 
are based on the photoevaporating flow models developed by [9, 4, 13, 16, 34]. 

The proplyds were discovered by Laques and Vidal in 1979 [17]. The term proplyd 
was coined by C.R. O’Dell to name the 6 nebulosities called LV objects identihed by these 
authors. The proplyds were identihed as young stars with circumstellar clouds photoionized 
from the exterior by 0^ Orionis C by means of VLA observations [6, 10]. They were imaged 
in 1993 with HST to show, for the hrst time, their morphology: silhouettes, bow shocks 
and even circumstellar disks [24, 25]. Indeed they were well understood to be a compact 
disks around YSOs with these excellent imagery. Recently the proplyds are designated 
acoording the notation of O’Dell & Wen [25]. This designation system is based on the 
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proplyds ecuatorial (J2000) coordinates (1” in declination and 1.5” in Right Ascension). By 
example, the proplyd 231-739 has coordinates (q:=5^ 35™ 23.123^, 5=-5° 27™ 38.9^). 

The proplyds have M of 10“® - 10“’^ yr“^ [6, 14], radial velocities of 24 to 30 km 
s“^ [14], diameters of 10“"^ - 10“^ pc [28] and disk masses of 0.005 - 0.02 [2, 16]. These 

values are estimated from extinction measures of the silhouette disks on the background 
HII region emission [36] or from CO milimetric observations of the thermal dust emission 
from their disks [2]. These latter observations give evidence of accretion disks surrounding 
young stellar objects ([34] and references therein). In the Orion Nebula (M 42, NGC 1976) 
there are ~ 150 proplyds [27], 21 present microjets (see Figure 7a of [3]) and 15 appear in 
pure silhoutte; YSOs and circumstellar disk (see Figure 7a of [3]). A recent review about 
the Orion nebula is presented by O’Dell [29]. 


2. ASTROBIOLOGICAL ASPECTS 

Given that in the protoplanetary disks, planet formation can occur and, as a 
consequence, it is possible that the formation of life can also occur there, the interest 
in astrobiological studies of the proplyds stems from the fact that proplyds unveil the 
otherwise ellusive protoplanetary disks. In this sense, proplyds are an important laboratory 
for determining, at high spatial resolution, several properties of the protoplanetary disks 
such as: masses, velocities and extensions. The proplyds are also useful in statistical 
studies on these issues. Furthermore, the protoplanetary disks revealed by the proplyds 
are subject to the strong UV flux of the exterior star that makes them visible. Thus, the 
photoevaporation of the protoplanetary disk of the proplyd could inhibit planet formation or 
change its conditions. Recent studies [2, 36] have shown that planet formation in irradiated 
protoplanetary disks as proplyds can only occur in special situations and conditions: disk 
mass >0.13 and dust (silicates-l-ices) particle radius > 5 fim. 



6 


Furthermore, Throop et al. [36] through Ha and Paa HST images, obtained an 
extintion curve of the biggest (in disk) proplyd 114-426 hnding that a theoretical extintion 
curve dominated by dust grains >5 /rm hxed well the observations. They also present 
1.3mm OVRO observations suggesting that the dust growth reaches radius values of few 
mm, and developed a numerical model to explain the behaivor and evolution of dust grains 
in proplyds irradiated externally including photo-destruction process. With this, they have 
shown that after 10^ yr, small grains are entrained in the photoevaporative flow resulting 
in disks with maximum sizes of 40 AU. In less than 10^ yr, in disk masses as large as 0.2 
the grain growth reaches a radius of meters at 10 AU of the YSO and of 1 mm at 500 
UA from it. These sizes allow the dust to resist the photoevaporation process. By 10® yr, 
nearly all ice and gas are removed by photosputtering inhibiting the Kuiper belts formation 
and leaving only the possibility of rocky planet formation. Therefore, it is possible that 
in the environments of star forming places like Orion, Jupiter-like planets are not formed 
in the standard way because it would require 10® to lO’^ yr ([33] and references therein). 
Such Jovian planets could be present if they form on 10^ yr time-scale [5] in disk masses 
> 0.13 Since most of the stars seem to be formed in large and dense clusters such 

as Orion, this leads to the conclusion that planet formation models should be revised in 
order to include the destructive effects of the UV flux of newly formed massive stars in star 
forming environments. The possibilities for planetary system to form in the Orion Nebula 
are described by Throop et al. [37]. 

In summary, in enviroments with external hot stars like the proplyds in the Orion 
Nebula: 1.- The formation of Jovian planets and kuiper belts objects is very dihcult, 
terrestrial planets remain unaffected and posible. 2.- The solar systems formed here must 
be different from our solar system, in fact, the standar solar nebula models can not be 
applied here. The important parameter to determinate the planet formation time scales is 
M. Indeed, the typical values suggest life time for the protoplanetary disks of 10® yrs. This 



last involve that in many of this proplyds, the planet formation can not occurr. Maybe the 
15 proplyds suggested by Bally and co-workers (see [3] and Bally in this symposium) are 
factible places for planet formation. 


3. FABRY-PEROT STUDY 

With Fabry-Perot (FP) interferometry and long-slit spectrographs, it is possible 
to study the kinematics of every astronomical object. Both techniques are useful and 
complementaries. The fundamental difference between both instruments is that with a FP 
is possible to study the whole region because the instrument cover all the held although 
the information is concentrated in a small part of the electromagnetic spectrum. In the 
spectrograph occurs the opposite: it is possible to obtain a spectrum but in a very limited 
zone of the object. Due to the fact that Orion is a nebula rich in interstellar phenomena 
and that it displays a gradient in the radial velocity across its area, we can, with FP 
techniques, identify the different velocity components of a radial velocity prohle, in this 
case, the HII region velocity component and the proplyd velocity component. Spectroscopic 
studies in [OIII] of the proplyds reported in [17] were made by [21, 19, 20, 22, 12] using slit 
spectrometers. Keck high resolution slit spectroscopy of proplyds at several ions including 
Ha, [SII] and [OIII] is presented by [14]. On the other hand, [26] present a Fabry-Perot 
kinematical study at [OIII](5007A) and [SII](6731A) of the inner region of the Orion 
Nebula, discovering new high velocity features and studying the redshifted and blueshifted 
emission of some proplyds. With these spectroscopical studies it is possible to determine 
the flow velocity, and from the [SII] lines and/or Ha emission measure, the density of the 
IF [1]. Knowing the systemic velocities and the radial velocity profiles of the proplyds 
we can compare the observations with the existing theoretical models in order to present 
arguments in favor of the existence of accretion disks in them, and know the behavior 



between the interaction of the photoevaporated material and the photoionizing winds. Also, 
from the determination of the velocity of the gas we can derive the M in combination with 
observations of the surface brightness in Ha [28]. 

This paper and [7] complements the previous works [31, 32], where a study of the 
kinematics of the gas, Herbig-Haro (HH) objects, and jets at large scales in the Orion 
Nebula are presented. These works between other things, find that HH 202, HH 203 and 
HH 204 are part of a large bipolar outflow centered near the E-W jet discovered by [26]. 
This last is suggested because HH 202 form part of a big blue-shifted lobe, whereas HH 
203-204 form part of a other big lobe that present red-shifted emission in our data cubes 
(e.g. see last panels on Figures 5 and 6 of [31]). Furthermore, recent proper motion study 
of these HH objects suggest that they emerge of the same source [8]. 


3.1. Observations, technique and results 

Fabry-Perot data cubes were obtained from November 30 to December 5, 1996 with the 
“PUMA” scanning Fabry-Perot spectrograph [30] at the f/7.5 Cassegrain focus of the 2.1m 
telescope of the Observatorio Astronomico Nacional at San Pedro Martir B.C., Mexico. See 
[31] for further details. The data reduction was made using the software CIGALE [18]. The 
results are presented in detail by [7]. 


3.1.1. Unsharp-masking technique and proplyds identification 

We were able to identify the proplyds: 82-336, 158-323, 158-326, 159-350, 163-317, 
167-317, 170-337, 177-341 and 244-440 from our raw FP data cubes (see Figure 1 of [31, 
32]. These proplyds show a conspicuous appearance in some of the velocity maps: they are 
detected as bright, point-like nebulosities. In order to improve the detection of proplyds 
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Table 1: Positions and velocity extent of the identifed proplyds 


Proplyd Name " 

x*- 


region 

(km s“^) 

^proplyd 

(km s~^) 

Vi„ 

(km s“^) 

FWHM proplyd 

(km s“^) 

FWHM corrected 

(km s“^) 

82-336 

193 

251 

13 

33 

37 

18.92 

6.5 

158-323 (LV 5) 

254 

228 

9 

34 

32 

18.92 

6.5 

158-326 (LV 6) 

255 

235 

11 

33 

- 

28.38 

20.14 

159-350 

251 

275 

9 

37 

33 

18.92 

6.5 

161-314 

246 

211 

9 

27 

14 

18.92 

6.5 

161-324 

248 

231 

9 

34 

- 

9.46 

17.62 

163-317 (LV 3) 

242 

218 

12 

24 

- 

18.92 

6.5 

166-316 

236 

212 

11 

27 

- 

18.92 

6.5 

167-317 (LV 2) 

231 

218 

12 

28 

- 

18.92 

6.5 

168-326 (LV 1) 

228 

235 

9 

36 

- 

28.38 

20.14 

170-337 

224 

253 

13 

22 

- 

28.38 

20.14 

176-325 

210 

231 

9 

36 

32 

18.92 

6.5 

177-341 

205 

259 

13 

32 

28 

37.84 

32.12 

180-331 

197 

241 

12 

38 

24 

18.92 

6.5 

197-427 

202 

295 

5 

30 

- 

18.92 

6.5 

244-440^= 

34 

361 

19= 

23= 

- 

47.30= 

42.86= 


"According with the notation of O’Dell Wen [25] 
^Pixel position in our data cubes 
"See the text 


from our velocity maps we applied an “unsharp-masking” process in order to isolate the 
small diameter emission structures from the bright diffuse HIT region (as it will be described 
below). This technique is also useful for detecting extended features embedded in bright 
emission (e.g. [23]). It consist in creating an out-of-focus image which is used as a spatial 
frecuency hlter, to enhance rendition of any hne details, while at the same time reducing 
gross density variations on the resulting hltered image. Indeed, this process improves the 
detection of proplyds. The proplyds: 161-314, 161-324, 166-316, 168-326, 176-325, 180-331 
and 197-427 became noticeable only after applying the unsharp-masking process while the 
proplyds already detected in the raw data cubes were more easily distinguished. 
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Figure 1 shows the resulting FP velocity map at V/^eiio = - 127 km s“^ (presented in 
[31, 32]) after applying the unsharp-masking process. In this map we have marked the 
position of the identihed proplyds that, as seen in this hgure, are quite conspicuous. The 
unsharp-masking process has been carried out in the following way: hrst, we have done a 
spatial, Gaussian smoothing (with a = 3 pixels or 1.77 ") that smoothes small diameter 
features (such as the proplyds and thin hlaments). Then, we subtracted the smothed 
velocity maps from the original velocity maps. The result is that the small diameter features 
are sorted out whereas the diffuse, extended emission (such as the foreground HII region) is 
subtracted. In that way we were able to identify the proplyds reported in [28]. The pixel 
coordinates of the identihed proplyds are quoted in Table 1 where the proplyds’ names 
according to O’Dell & Wen are given in column 1, whereas the proplyds pixel coordinates, 
X and y, are given in columns 2 and 3 respectly. 


3.1.2. Proplyds’ radial velocity profiles and profile decomposition 

Once we have obtained the pixel coordinates of the identihed proplyds, we proceeded 
to extract radial velocity prohles integrated over boxes of 2 x 2 pixels centered at the 
pixel positions reported in Table 1. For the obtainment of the radial velocity prohles we 
use only the original Ha data cubes (i.e., before applying the unsharp masking process). 
The radial velocity prohles obtained in this way are contaminated by the bright HII region 
velocity component (proplyd-|-HH region prohles). We preferred to subtract the HII region 
contribution point by point taking advantage of the fact that our FP data cubes give also 
the velocity prohles of neighboring regions outside the proplyd position. In this way, in 
order to know the contribution of the HII region emission, we extracted radial velocity 
prohles of several zones close to the proplyd location. This allows us to identify the average 
intensity and velocity of the HII region velocity component (HII region velocity prohle). 
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Finally, using this information, we fitted in the proplyd+HII region prohle, the HII region 
velocity prohle and a remaining velocity component that we identify as the proplyd velocity 
prohle. The main problem and difficulty in the kinematical studies of the proplyds resides 
in separating the velocity component of the proplyd from those of the HII region because 
the intensity of the HII region with respect to the proplyd is higher. Furthermore, both the 
proplyd and the HII region can have similar or same radial velocity making the task quite 
difficult. However, the use of FP techniques that cover all the held around the proplyd, 
allow us to get the radial velocity prohle of the HII region quite accurately. Consequently, 
we were able to obtain the proplyd prohles with an accuracy of ± 5 km s“^ in the peak 
velocities and ± 10 km s“^ in the FWHM’s. 

Figures 2, and 3 show the radial velocity prohles obtained at the identihed proplyds 
positions, as described above. The radial velocity prohles shown in these hgures correspond 
to the proplyd+HII region prohles. The velocity components of the HII region and those 
assumed to be of the proplyds are also shown in this hgure. Table 1 also quotes the 
heliocentric peak radial velocities of the contaminating HII region and the proplyds, as well 
as the FWHM of the proplyds (columns 4, 5 and 7 respectively). The HII and proplyd 
velocity components are in agreement with typical values cited in the literature (e.g. [28, 
14]). As it is discussed in the study by [26], many proplyds present redshifted high velocity 
hows. For 159-359, 82-336, 180-331, 175-325 and 161-314 we also found weak blueshifted 
components (column 6 in Table 1) probably related to the HII region emission (Figures 2 
and 3). About these components, a better spatial resolution studies are needed to clarify 
their origin. Unfortunately, we can only compare our results with the velocity prohles 
obtained from spectroscopic Ha data of 3 proplyds (170-337, 177-341 and 244-440) of [14]. 
We hnd a good agreement with 170-337 for the proplyd peak velocity (22 vs 20 km s“^) 
while for 244-440 and 177-341 we hnd that the proplyd peak velocity is diherent from the 
one reported by [14] of 23 vs 10 km s“^ and 32 vs 22 km s“^, respectively, in spite that 
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those proplyds appear bright in our data cubes. The values reported in Table 1 for 244-440 
cover all its area. Indeed, a more detailed inspection (see [7]) suggest that the center (1x1 
pixels) of the proplyd have Vheiio = 16 km s“^, whereas the regions around it have values 
of 18-23 km s“h The detection of proplyds in our [Nil] data cubes is poor. Indeed, we 
only identify two of them; 170-337 and 244-440. Howewer, they have radial velocities of 29 
and 12 km s“^ that are in agreement with the results reported by [14] of 33 and 10 km s“^ 
respective! ly. It is very important to carry out 3-dimensional spectroscopy especially in 
Ha, [01] and [Nil], with better spectral resolution but covering the entire held of view, in 
order to verify whether there is agreement or not with the results reported here. 


Our results do not allow us to study the proplyd prohle, including their diferent 
components and even their forms (not necessarily Gaussian). We can only hud other 
velocity components different from the HII region, proposing that these velocities are 
related with the proplyds. Table 1 lists the proplyd peak velocities as well as the velocities 
of the contaminating HII regions obtained from our FP velocity prohles. As one can see 
from this table, all the proplyd peak velocities are redshifted relative to the HII region 
indicating that all the proplyds studied in this work are located in the same gaseous sheet 
behind the HII region. In addition. Table 1 reports the proplyd velocity widths uncorrected 
and once corrected from the instrumental and thermal widths. These widths have values 
of 10 and 20 km s“^ respectly, similar to the presented by [14]. This table also quotes the 
Ha velocity widths (FWHM) obtained in this work, corrected for the instrumental function 
and thermal broadening (column 8). An inspection of the values of this quantity for all the 
proplyds we study here shows that the velocity widths are either smaller or comparable 
to the ones reported by [14], and that the Ha velocity widths do not seem to vary as a 
function of distance to the exterior ionizing source. Consequently, we expect that the initial 
flow velocity is comparable to the one found by these authors, i.e. vo=13 km s“^. With 
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this value and taking ng, tq from the literature [18] (as is discussed and presented in de la 
Fuente et ah [7]), we found M values between 10“® to 10“’^ yr and r ~ 10"^-10® yr. 
By example, the proplyd 161-314 have M ~ 2x10“^ yr and r ~ 5x10® yrs. Thus, 
our kinematic work conhrms the conclusions on time scales for disk destruction of 10® yr 
[13]. Futhermore, our r values suggest that in our proplyds sample, the planet formation 
could be difficult and indeed, implies the necessity to revise the models of planet formation 
as discussed by [36, 37]. It should be quite valuable to obtain FP observations at [01](6300 
A) in order to derive with better accuracy, the M of each proplyd. Futhermore there is 
some evidence th! at vq is traced by this ion [14]. 


MR wishes to acknowledge the hnancial support from DGAPA-UNAM via the grant 
IN104696. E de la F wishes to acknowledge hnancial support from CONACYT-Mexico 
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Fig. 1.— Unsharp masking Ho; velocity map at heliocentric velocity = -127 km s“^. The 
original image published by Rosado et ah [31, 32] is shown in the insert. The identihed 
studied objects are marked. 

Fig. 2.— Radial velocity prohles of the detected proplyds 82-336, 158-326, 159-350, 161-314, 
161-324 and 163-317 (longest component). The middle size component stands for the HII 
region. The smallest one represents the suggested proplyd. Some prohles feature an extra 
blueshifted component. The squares represent the ht of these components. 

Fig. 3.— Same as Figure 2 for proplyds 166-316, 167-317, 168-326, 170-337, 176-325, 180-331, 
158-323, 177-341 and 197-427. 



